Nuclear quadrupole resonances in compact vapor cells: 
the crossover between the NMR and the nuclear quadrupole resonance interaction regimes 
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We present an experimental study that maps the transformation of nuclear quadrupole resonances from the 
pure nuclear quadrupole regime to the quadrupole-perturbed Zeeman regime. The transformation presents an 
interesting quantum-mechanical problem, since the quantization axis changes from being aligned along the axis 
of the electric-field gradient tensor to being aligned along the magnetic field. The large nuclear quadrupole shifts 
present in our system enable us to study this regime with relatively high resolution. We achieve large nuclear 
quadrupole shifts for 1 — 3/2 131 Xe by using a cube-shaped 1 mm 3 vapor cell with walls of different materials. 
The enhancement of the NQR shift from the cell wall materials is a new observation that opens up an additional 
adjustable parameter to tune and enhance the nuclear quadrupole interactions in vapor cells. As a confirmation 
that the interesting and complex spectra that we observe are indeed expected, we compare our data to numerical 
calculations and find excellent agreement. 

PACS numbers: 32.60.+i, 33.25.+k, 76.60.Gv 



I. INTRODUCTION 

Any atom that has a nuclear spin I > 1 has a nuclear elec- 
tric quadrupole moment, whose interactions with electric field 
gradients can cause shifts of the nuclear magnetic energy lev- 
els. There is a large body of literature on the interactions of 
nuclear quadrupole moments with electric field gradients. As 
far back as the 1950s, studies were performed in crystals both 
in the regime where the nuclear quadrupole interaction caused 
weak perturbations to the nuclear magnetic resonance (NMR) 
spectra [1], as well as in the pure nuclear quadrupole reso- 
nance (NQR) regime, where little or no Zeeman interaction 
was present [2]. Solutions for the transition energies between 
nuclear spin sublevels were found for both regimes by use 
of perturbation theory (for a review, see y[]) by aligning the 
quantization axis along the principal axis of the electric field 
gradient tensor in the NQR regime and along the axis of the 
magnetic field in the NMR regime. As with these first experi- 
ments, most NQR studies have either been distinctly in either 
the NQR or the NMR regimes. To our knowledge, prior to 
our work reported here, the transformation of the NQR spec- 
tra from the NQR to the NMR regimes had not been observed 
experimentally. 

Cohen-Tannoudji first suggested that in addition to arising 
from electric field gradients from ionic bonds in a crystal, 
quadrupolar coupling could occur between nuclei and electri- 
cal field gradients present at the nucleus during wall collisions 
for atoms in vapor cells Q4|] . Since then, nuclear quadrupole 
resonances in vapor cells have been studied for many systems 
including I = 3/2 2O1 HgJlJE0JEJl, I = 9/2 83 Kr d, 
I = 3/2 131 Xe liltilllltllllltll, and / = 3/2 21 Ne 
Inll . Much of this work has been of basic interest from a 
fundamental physics standpoint JH 0, 0, US [Tl ljl^. [l3i [Til 
[l5l,[l6ll and for tests of fundamental s vmmetrie^llSlM fl7l 1 1 8ll . 
There have also been proposals for using these systems for the 
practical application of rotation sensing II 191, 12011 . Changes to 
the NQR shifts in the crossover regime could lead to system- 
atic errors in precision measurements and offsets in rotation 



sensors. 

For much of the NQR work that has been performed in va- 
por cells, the NQR lines were not clearly resolved and the 
NQR splitting caused a slow beatingor a nonexponential de- 
cay of the nuclear polarization JH |a OH EH- Th e beat 
frequencies depended on the orientation of the cell symmetry 
axis in the magnetic field, and went to zero at the magic an- 
gle of 54.7 J5|], which indicated that the axis of symmetry of 
the electric field gradient was aligned along the cell symmetry 
axis. 

In a series of papers from Happer's group at Princeton 
lfl2[ [l~3l Wu et al. saw much stronger interactions such 
that the NQR lines could be clearly resolved by using highly 
asymmetric cells. They performed a detailed perturbation- 
theory solution for the NQR shift in the NMR regime, ac- 
counting for pressure-dependent diffusion and cell shape, and 
formulated the results to give a microscopic description of 
the interaction lfl3ll . Ignoring complications from diffusion, 
they expressed the NQR shifts for the [ — 3/2) < — 1/2 1 and 
1 1/2} (3/2 1 coherences as 
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which is an integration of the nuclear quadrupole interaction 
over the cell walls. Here v is the atom velocity, (9) is the 
mean twist angle per wall adhesion, S is the cell surface area, 
V is the cell volume, / is the nuclear spin, and ijj is the an- 
gle between the local surface normal (directed out of the cell) 
and the magnetic field. Here we put (9) within the integral 
to allow for the possibility that the cell walls are of differ- 
ent materials. Integrating Eq. Q]for a cylindrical cell gives 
AS1 = ±AJ7oP2 (coscp), where <p is the angle between the 
cell symmetry axis and the direction of the magnetic field 
andP 2 (x) = i(3x 2 -l). AO = vA (9) S/4V is pro- 
portional to the atom velocity, the surface to volume ratio of 
the cell, and an asymmetry parameter, A, which goes to zero 
when the cell height and diameter are equal. Wu et al. ver- 
ified their theory with detailed experiments and determined 
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(9) = 38(4) x lCT 6 rad for pyrex f\M . Experiments later 
performed by Butscher et al. lfl5ll revealed similar behavior. 

All of the experiments described so far were performed at 
magnetic fields high enough that the nuclear quadrupole in- 
teraction was well described as a perturbation to the magnetic 
Larmor resonances. Appelt et al. performed studies on 131 Xe 
in the limit of zero applied magnetic field and measured devi- 
ations from Berry's adiabatic phase under rotation [16]. They 
solved the 1 = 3/2 Hamiltonian by including terms for the 
NQR interaction and spatial rotations and showed that mix- 
ing of the nuclear spin sublevels through rotation makes all 
six transitions between nuclear sublevels allowed (Am =1,2, 
and 3). 

All of this work was either distinctly in the NMR JT^, [H 
[H [H [H QJ] or the NQR d] regime. Here we bridged 
the two regimes by continuously tracing the transformation 
from the pure NQR regime to the quadrupole-perturbed Zee- 
man regime. We achieved a large NQR splitting not by using 
geometrically asymmetric cells as Wu et al. 01211 did , but 
by using a 1 mm 3 cubic cell with walls of different materi- 
als. A cubic charge distribution would not ordinarily cause 
an NQR splitting [7], but the microscopic surface interactions 
with the different wall materials lower the symmetry of the 
system. The small cell size enhances the NQR splitting be- 
cause the shifts are proportional to the surface to volume ratio. 
The enhancement of the NQR shift from the cell wall materi- 
als is a new observation that gives one an additional adjustable 
parameter to tune the electric field gradient and the NQR in- 
teraction in vapor cells. As a confirmation that the complex 
spectra that we observe are expected, we compare the trans- 
formation of the resonance lines to theoretical models and find 
excellent agreement. 



II. EXPERIMENT 
A. Techniques and Apparatus 

Fig. [TJ\ is a schematic drawing of our apparatus. The mi- 
crofabricated sample cell of volume 1 mm 3 was etched in sili- 
con and sealed with pyrex [21]. The cell contains 87 Rb, buffer 
gases of N2 and Ne at 10 and 600 torr, respectively, and 10 torr 
of Xe gas at natural abundance. Xe has two active NMR iso- 
topes: spin- 1/2 129 Xe (26.4% abundance) and spin-3/2 131 Xe 
(21.2% abundance). The cell is cubic and has four silicon 
walls and two pyrex windows. The cell is heated to 145 °C 
and mounted at the center of a set of three orthogonal mag- 
netic coils. The coils are surrounded by a four-layer magnetic 
shield (one layer shown in Fig. UK). A circularly polarized 
laser beam optically pumps and probes the Rb atoms through 
the pyrex cell windows along the z-axis. The Rb polarizes the 
Xe atoms through spin-exchange optical pumping Q. The 
Rb also functions as a magnetometer and is used to sense the 
magnetic fields generated by the Xe atoms ll23l l24ll . 

We use a field switch technique to initiate precession of the 
Xe atoms. During the pump phase, the Xe polarization builds 
and reaches steady state. At the start of the probe phase, a DC 
magnetic field, Bq, is turned on in the y—z plane and the Xe 
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FIG. 1: (A). The apparatus. The 1 mm cell is roughly centered 
on an orthogonal three-axis set of magnetic coils. A laser beam of 
maximum power 3mW enters the cell. The transmitted power is 
detected by a photodiode. (B). The coordinate system: the z-axis 
coincides with the light propagation direction, k. The magnetic field 
direction is in the y-z plane rotated by an angle tp from the i-axis. 



atoms start to precess. The angle of Bq with respect to the 
z-axis (the cell symmetry axis), ip, is varied, depending on the 
experiment. An AC magnetic field, Bac, of rms amplitude 
~ 1//T and frequency ~ 2 kHz drives the Rb atoms and is ap- 
plied along the x-axis. This AC drive also references a lock- 
in amplifier that measures the modulation of the transmitted 
power at the Rb drive frequency. The applied field geometry 
is shown in Fig. [Tj3. We are also able to observe signals in 
many other field configurations, but we have found this con- 
figuration to give the best signal-to-noise ratio. Our geometry 
for pumping and probing is very similar to the technique used 
by Volketal. Q. 

After a field switch, a free induction decay (FID) signal is 
observed at the output of the lock-in amplifier. Fig. [2] shows 
an FID signal and its Fourier transform. In most cases, we 
used an acquisition time of 32.8 s on our spectrum analyzer, 
which gave a frequency resolution of 30.5 mHz. Since the 
acquisition time was much longer than our typical T2 time 
of ~ 5 s, we compromised on signal-to-noise ratio to achieve 
higher frequency resolution. 

One factor that complicates our estimating the field mag- 
nitude and angle is the relatively large field generated by the 
Rb atoms as sensed by the Xe atoms B26I1 . the magnitude of 
which is Bfib = M-B^pf^flfc-P- Ms is the Bohr magneton, 
k = 730 is the hyperfine contact enhancement factor Il27ll . /zo 
is the magnetic constant, nut, is the Rb density, and P is the 
Rb polarization. At our laser intensity of 300 mW/cm 2 and 
temperature of 145 °C, we measure But = 200 nT, which 
corresponds to P = 30 %. 

To simplify controlling the total field in the presence of this 
large offset field, for most of our measurements we divide our 
applied field into two components such that our total field is 
B to t = Bo + Bin, + B comp . We set the compensation field, 
B comp , equal to —Bub, such that we can determine our total 
field as sensed by Xe from the variable component Bo alone. 
Offsetting Bin, is made easier by having a high pumping rate 
so that Bri> II k is independent of the angle of B Q . Then B^ 
can be nulled out with a constant field parallel to the direc- 
tion of light propagation. This approach simplifies controlling 
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both <p and the total field as the angle or magnitude of Bq is 
varied. 



B. Measurements 

We concentrate our measurements on two things: first, the 
131 Xe NQR shift Ail versus angle p, and second, the transfor- 
mation of the energy shifts as Bq is swept from zero through 
the NQR-dominated regime and into the NMR-dominated 
regime. To measure A57 (Eq. 1) versus ip, Bq was kept near 
0.8 fiT. At high enough field, the NQR splitting depends on 
the field angle but not on the field magnitude. For this mea- 
surement, we did not apply a compensation field as described 
above, but rather included B^h in our calculation of (p. The 
data were consistent with Bri, \\ k. We measured the 129 Xe 
and 131 Xe spectra versus angle at two different laser powers: 
0.6 mW and 3 mW (60 and 300 mW/cm 2 ). We determined fre- 
quency differences between the outer two 131 Xe resonances 
using curve fitting and took Af2 to be half of the difference. 
The inset in Fig. [2j3 is a plot of Afi versus cos p for two dif- 
ferent laser powers. Only field angle values where the triplet 
is resolved for curve fitting are included. The parabolic fit 
of the data agrees well with Aft — at the magic angle of 
p = 54.7° (cos ip = 0.578). 

To make quantitative estimates of the mean rotation angle 
per wall collision, we apply Eq. Q]to our cubic cell. Assuming 
(9) for pyrex can be expressed in terms of (9) for silicon as 
(9 p ) = (9 S ) +5 one finds A51 = Af2oP2 (cos p), where Ai7o 
reduces to vd/L. L is the length of a side of the cube, and v is 
the mean thermal velocity. Fitting this expression to the data 
in the inset of Fig. G£, we find Afi = 2tt x 0.39 Hz. With 
L = 1 mm, and v = 281 m/s, we find i5 = 8.7 /xrad. Using Wu 
et al.'s measurement of (9 p ) — 38 /irad for pyrex 01411 . we find 
(9 S ) = 29 /irad for silicon. 

To measure the spectra versus field magnitude for a fixed 
field angle, we carefully offset Bm, by applying a compensa- 
tion field along k. We collected multiple spectra versus field 
magnitude at several field angles. Data are presented in Fig. 
|3]for <p =22° and 39°. The plots are three-dimensional - the 
vertical axis is the measured frequency, the horizontal axis is 
the applied magnetic field (not counting B comp ), and the sym- 
bol size is proportional to the signal amplitude. Each vertical 
line represents an individual frequency spectrum collected at 
a fixed magnetic field. The black lines show the transition 
frequencies for 129 Xe versus magnetic field. Since 129 Xe has 
no nuclear electric quadrupole moment, the transition energy 
is linear in applied field. The gray lines represent the energy 
differences between the four nuclear sublevels of 131 Xe found 
numerically, as discussed below. 



III. CALCULATIONS 

In the limit that the NQR shift is either much smaller or 
much larger than the Larmor frequency, perturbation theory 
solutions accurately predict the transition frequencies HI [3 • 
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FIG. 2: Examples of experimental data. (A). A free-induction decay. 
The signal is proportional to the amplitude of the modulation of the 
photodiode signal caused by the Rb precession and is presented in 
volts. The magnetic field during the probe phase was 275 nT, and ip 
was ~ 70°. Beating from the 131 Xe triplet is seen in inset, which 
shows the residuals from a fit of a damped sine wave to the data. (B). 
The Fourier transform of the data in (A). Here the signal is presented 
as a power spectral density in the units of root-mean-squared volts 
in a 1 Hz bandwidth. The inset is a plot of the NQR splitting versus 
cos (p. The solid symbols are data collected at 300 mW/cm 2 , and the 
open symbols are data collected at 60 mW/cm 2 . For the low (high) 
power data, the Rb field was 0.12 /iT (0.19 fiT). For comparison, the 
applied field was ~ 0.8/iT. There is more distortion in the parabolic 
shape of the curve for the low power data - particularly when the 
applied field was orthogonal to the light propagation direction where 
the approximation that But || k would not hold as well. The dashed 
line is a fit to the high-power data, which gives a value of A£lo/2n = 
0.39(1) Hz. Our value for Afi is ~ 43% larger than the largest 
value reported by Wu et al. fbUl . 



When the NMR and NQR interactions are of comparable size, 
a more involved solution is required. It is possible to solve the 
system analytically by diagonalizing the full Hamiltonian to 
find the transition energies as well as the transition amplitudes 
ll28ll . but it is involved - particularly given the dynamics and 
the angular sensitivity of the Rb magnetometer. Finding the 
full analytical solution is made dramatically more accessible 
by taking as a starting point a general solution that has already 
been developed and is available online [29]. A full calcula- 
tion yielding predictions for the line amplitudes that will give 
more insight into the physics is underway and will be reported 
elsewhere. For this work, we compare our data to numerical 
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FIG. 3: Spectra versus applied field for p = 22°(A) and tp — 
39°(B). The grey lines represent the six frequency differences be- 
tween the four nuclear sublevels computed numerically (see text). 
The black lines mark the 129 Xe transition versus field. Each ver- 
tical line represents a different frequency spectrum acquired at the 
specific magnetic field that it intersects on the horizontal axis. The 
symbol size is proportional to the signal amplitude, and the same 
amplitude scale was used for both angles of p. The laser intensity 
was 300 mW/cm 2 and B Rb was 0.19 fiT. The inset in (A) shows the 
131 Xe nuclear energy levels in the low-field NQR regime. Since the 
energies of the nuclear sublevels go as m 2 at zero magnetic field [1], 
the 1+3/2} and | — 3/2) states are degenerate at zero field, as are the 
| + l/2) and | — 1/2) states. The inset in (B) shows the nuclear energy 
levels in the high-field NMR regime. 



quencies between the four nuclear states. At high magnetic 
fields in the NMR regime, there are three lines with equal 
slope, corresponding to Am = 1 transitions (|— 3/2)(— 1/2|, 
| - 1/2) (1/2 1, and 1 1/2) (3/2 1), two lines with double the slope 
for Am = 2 transitions (|-3/2)(l/2| and |-l/2)(3/2|), and 
one line for the Am = 3 transition (|-3/2)(3/2|). Note that 
these transition frequencies are plotted with no insight into the 
transition amplitudes, and in fact, only the Am = 1 transi- 
tions are allowed at higher fields. At very low magnetic field, 
Am = 1,2, and 3 transitions are also seen, but they do not 
correspond to the same Am = 1,2, and 3 transitions seen at 
higher field because the quantization axis rotates as a function 
of magnetic field thus transforming the states. 

The observed 131 Xe lines agree well with the predicted fre- 
quencies even as the magnetic field goes to zero. At low field, 
the transitions are unresolved, and in most cases it is difficult 
to assign features to Am = 2 or Am = 3 transitions. Per- 
haps the best resolved spectrum at low field, where we would 
expect strong mixing of the lines, was collected for 28 nT and 
ip = 39°, which corresponds to the second spectrum from the 
left in Fig. [3jJ. This spectrum is shown in Fig. [4] The pre- 
dicted transition frequencies are marked. One of the Am = 2 
transitions is clearly visible and is one-tenth as strong as the 
Am = 1 transitions. The other Am = 2 transition and the 
Am = 3 transition are not visible. We cannot put limits on 
the strength of the Am = 3 transition, since the frequency 
where it would appear is obscured by the wings of other, much 
stronger transitions. 



IV. DISCUSSION AND OUTLOOK 

One point requiring further investigation relates to the am- 
plitudes of the 131 Xe lines. Whereas the 129 Xe line ampli- 
tudes vary by about 10% as Bq is varied, the 131 Xe line am- 
plitudes vary by much more - sometimes jumping up by a 



calculations using a Liouvillian approach described in detail 
by Bain 13011 . The method is relatively straightforward to use 
for calculating the transition frequencies for an arbitrary spin 
nucleus in an arbitrary magnetic field and electric field gradi- 
ent, but it does not predict the transition amplitudes and it is 
also not very transparent. The details of the calculation are 
outside of the scope of this article, and we refer the interested 
reader to the original work [30], which gives the recipe for the 
calculations. 

Three parameters enter the calculation: the angle ip, Afio, 
and an asymmetry parameter, 77, which is zero in the case of a 
cylindrically symmetric electric field gradient. We set r\ = 
for our simulations and for Af2o we use our measurement 
from the inset of Fig. |2jJ. We used our estimates of ip from our 
coil calibrations and cell orientation assuming the symmetry 
axis of the electric field gradient to be along the cell symme- 
try axis. We conservatively estimate an upper limit of 5° for 
angular misalignment of the cell axis from the magnetic field 
axis. 

The gray lines presented in Fig. [3]are the six transition fre- 
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FIG. 4: A single spectrum collected with 28 nT and tp — 39°, which 
corresponds to the second spectrum from the left in Fig. [3] The pre- 
dicted transition frequencies are also shown for reference. The signal 
has the units of root-mean-squared amplitude for a 32.768 s averag- 
ing time. Given that the T2 time is about 5 seconds, we compromised 
on signal amplitude to achieve improved frequency resolution. 



5 



factor of 2 to 3 when the lines cross, and sometimes fading 
away and disappearing at a different field. The variations in 
line amplitudes for 131 Xe will be the subject of further study 
in conjuction with performing a full analytical solution. 

Our NQR shifts are large enough that we are able to mea- 
sure the decay rates for the individual lines. Like the line am- 
plitudes, the decay rates for the 131 Xe lines vary widely - es- 
pecially at the line crossings. The ratio of the decay rates for 
the |-3/2}(-l/2| and and |l/2)(3/2| coherences relative to 
the | — 1/2) (1/2 coherence does not agree with the 3:2 result 
predicted by Wu et al. for the NMR regime lfl3tl . The variation 
in the decay rates will also be an area of further study. 
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